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Abstract:  Statistical  and  spatial  analyses  were 
performed  on  existing  thaw-depth  data  sets  from 
a  variety  of  climatic  and  terrain  conditions  to  deter¬ 
mine  the  variability  of  thaw.  These  data  include 
permafrost  sites  in  Alaska,  temperate  sites  in  the 
continental  United  States  in  Minnesota  and  Wis¬ 
consin,  and  a  large-scale,  controlled  test  section 
with  no  solar  input  in  the  Frost  Effects  Research 


Facility  (FERF)  at  the  Cold  Regions  Research  and 
Engineering  Laboratory  (CRREL)  in  Hanover,  New 
Hampshire.  Results  can  be  used  to  spatially  dis¬ 
tribute  soil  properties  based  on  point  data  or  one¬ 
dimensional  models,  or  to  populate  sparse  data 
sets  with  terrain  properties.  They  are  also  useful  for 
analysis  of  impact  of  thaw  distribution  on  predic¬ 
tive  models,  such  as  for  predicting  vehicle  mobility. 
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Spatial  Analysis  of  Thaw  Depth 

ROSA  T.  AFFLECK  AND  SALLY  A.  SHOOP 


INTRODUCTION 

The  objective  of  this  work  was  to  develop  an 
approach  to  characterizing  the  spatial  variability  of 
thawing  soil.  Results  can  then  be  applied  to  spatially 
distribute  soil  properties  based  on  point  data  or  one¬ 
dimensional  models,  or  to  fill  sparse  data  sets  with  ter¬ 
rain  properties.  Results  are  also  useful  for  probabilistic 
analysis  of  the  impact  of  the  input  distribution  on  the 
confidence  of  predictive  models,  such  as  for  predicting 
vehicle  mobility  or  surface  runoff.  Thaw  depth  was  used 
as  the  variable  of  interest  because  of  its  clear  impact  on 
terrain  strength  and  hydrology.  Other  soil  variables  can 
be  analyzed  similarly.  A  comparison  among  these  data 
sets  also  indicates  the  impact  of  other  terrain  parame¬ 
ters,  such  as  vegetation  or  radiation,  on  the  thaw  varia¬ 
bility  for  a  particular  terrain  unit. 

A  search  and  review  of  existing  data  sets  having 
spatial  thaw  measurements  yielded  four  major  data  sets 
suitable  for  analysis.  These  data  cover  a  variety  of  cli¬ 
matic  and  terrain  conditions,  including  permafrost  sites 
in  Alaska,  temperate  sites  in  Minnesota  and  Wiscon¬ 
sin,  and  a  large-scale,  controlled  test  section  with  no 
solar  input  in  the  Frost  Effects  Research  Facility  (FERF) 
at  the  Cold  Regions  Research  and  Engineering  Labora¬ 
tory  (CRREL)  in  Hanover,  New  Hampshire.  Each  of 
these  data  sets  is  described  in  detail  below. 

Thaw  is  a  complicated  process,  influenced  by  the 
surface  energy  balance,  terrain,  and  soil  properties  at 
each  point.  Thaw  depth  is  the  thickness  of  the  thawed 
layer  from  the  surface  to  an  underlying  frozen  layer, 
which  typically  occurs  when  temperature  is  above  freez¬ 
ing  in  spring.  Terrain  with  a  thawed  layer  can  be  very 
saturated  due  to  snowmelt,  rain,  and  water  drawn  to 
the  soil  during  the  freezing  process.  Slopes  (aspect  and 
degree),  types  of  vegetation,  and  soil  properties  affect 
the  variation  on  thaw  depth.  With  all  these  variable  con¬ 


ditions,  one  would  expect  a  very  high  degree  of  hetero¬ 
geneity.  The  basic  question  is  to  define  the  spatial  vari¬ 
ability  of  depth  of  thaw  and  analyze  the  site’s  specific 
effects  on  this  variability.  The  Alaska  North  Slope  data 
sets  are  from  areas  of  continuous  permafrost  where  the 
underlying  frozen  layer  will  remain  year-round.  The 
bulk  of  the  measurements  were  taken  at  the  end  of  the 
summer  when  the  thaw  had  penetrated  to  the  maximum 
depth.  The  other  sites  are  located  in  areas  without  per¬ 
mafrost  where  all  frost  will  eventually  be  removed  as 
the  bottom  of  the  frozen  layer  meets  with  a  thawed  layer 
by  late  spring  or  early  summer.  The  data  sets  were  col¬ 
lected  over  a  wide  range  of  spacing  and  scales,  all  of 
which  affect  the  resulting  spatial  analysis.  This  is  an 
initial  effort  to  develop  variation  of  thaw  depth  for  a 
variety  of  conditions,  but  this  work  excludes  the  im¬ 
pact  of  surface  energy  balance  and  climatic  conditions 
on  the  thaw  variability. 

DATA  SETS 

Alaska  North  Slope 

Thaw  depth  (active-layer  thickness)  data  on  Alaska 
North  Slope  were  collected  for  Arctic  climate-change 
research  headed  by  Dr.  Frederick  E.  Nelson  for  Arctic 
System  Science/Land-Atmosphere-Ice  Interactions 
(ARCSS/LAII)  sponsored  by  the  National  Science 
Foundation  (Nelson  et  al.  1997).  The  National  Snow 
and  Ice  Data  Center  (NSIDC)  managed  and  archived 
these  data  sets. 

Permafrost  underlies  the  area.  Nelson  et  al.  (1997) 
used  Geographic  Information  Systems  (GIS)  to  map 
and  develop  an  integrated  response  of  thaw  depth  us¬ 
ing  air  temperature,  vegetation,  soils  and  soil  moisture, 
and  solar  radiation  for  these  sites.  Figure  1  shows  the 
locations  of  the  seven  sites  where  thaw  depth  was  meas- 
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Table  1.  Locations  of  the  seven  grid  sites  on  Alaska  North  Slope. 

Universal  Transverse 
Mercator  (UTM) 


Site 

Terrain 

description 

Coordinates* 

Zone  6t 

Latitude 

Lonqitude 

Northing 

Eastina 

Atquasuk 

Coastal 

70°27'N 

1 57°24'W 

7816500 

559200 

Barrow 

Coastal 

71°19'N 

156°36'W 

7913100 

585700 

Betty  Pingo 

Coastal 

70°17'N 

148°52'W 

7797700 

428900 

Happy  Valley 

Foothill 

69°06'N 

148°30'W 

7671500 

426000 

Imnaviat  Creek 

Foothill 

68°30'N 

149°30'W 

7612300 

405600 

Toolik  Lake 

Foothill 

68°37'N 

149°36'W 

7613900 

393500 

West  Dock 

Coastal 

70°22'N 

148°34'W 

7807300 

440800 

*  General  location  of  the  seven  sites. 

t  Universal  Transverse  Mercator  (UTM)  coordinates  for  the  southwest  corner 
(e.g.,  point  11  on  Fig.  2)  of  each  grid. 


ured  in  July  and  August  1995.  The  coordinates  of  the 
sites  are  given  in  Table  1 . 

Thaw  depth  was  collected  at  100-m  intervals  on  1-  by 
1-km  grid  (Fig.  2).  Thaw  depths  were  not  measured 
where  grid  points  intersected  lakes,  rocks,  or  roads. 
Thaw  depth  was  measured  using  a  rigid  metal  rod 
pushed  into  the  ground  to  refusal. 

Atquasuk  soils  are  mostly  sand  to  loamy  sand.  The 
soil  landforms  in  Atquasuk  are  grouped  into  three  cate¬ 
gories:  moderately  well  drained,  where  the  soils  have 
coarse  texture  and  low  concentration  of  organic  matter; 
somewhat  poorly  drained,  where  the  water  table  is 
below  the  surface;  and  very  poorly  drained,  where  the 
soils  have  considerable  amounts  of  organic  matter  and 
the  water  table  is  at  or  above  the  surface  (Everett  1 980). 
Everett  classified  the  soils  as  Pergelic  Cryaquepts 
(covers  most  of  the  area),  Pergelic  Cryohemists  and 
Cryosapriest,  Pergelic  Cryopsamments,  and  Pergelic 
Cryo fluents.  The  Atquasuk  area  is  considered  fairly  flat 


meadow  terrain.  Tussock  tundra  is  the  most  common 
type  of  vegetation  in  the  area.  A  detailed  description  of 
the  vegetation  near  Atquasuk  is  provided  in  Komarkova 
and  Webber  (1980). 

Barrow,  the  northernmost  point  in  Alaska,  is  domi¬ 
nated  by  wet  and  acid  soils.  As  described  in  Brown  et 
al.  (1980),  the  soils  in  the  area  are  formed  on  flat  to 
gently  sloping  topography  under  cold  and  moist  condi¬ 
tions  and  consist  of  a  considerable  amount  of  organic 
matter.  These  soils  have  a  three-layer  morphology:  an 
organic-rich  surface  horizon;  a  horizon  of  silty  clay  to 
silty  loam  textured  mineral  soil  with  variable  amounts 
of  organics;  and  an  underlying  horizon  of  a  perennially 
frozen  organic-rich  layer.  The  vegetation  of  the  gener¬ 
al  area  in  Barrow  consists  of  sedge  meadows  and  mosses 
(Brown  et  al.  1980). 

Betty  Pingo  is  located  in  a  wetland  terrain  with  min¬ 
imal  topographic  relief  (Hinzman  et  al.  1998).  The  soil 
in  Betty  Pingo  consists  of  silt  loam  soil  mixed  with  a 
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Figure  2.  Sampling  grid  for  the  seven  Alaska  sites. 


substantial  amount  of  organic  materials,  and  is  classi¬ 
fied  as  Euic  Pergelic  Cryosapriest  (Ping  et  al.  1994).  The 
vegetation  includes  sedge,  shrub,  lichens,  and  mosses. 

Happy  Valley  is  on  a  slight  slope  with  a  western 
aspect.  The  soil  is  classified  as  a  fine-loamy,  mixed 
Pergelic  Cryaquept,  which  is  composed  of  silt  loam  with 
organic  materials  and  a  layer  of  peat  on  top  (Ping  et  al. 
1 994).  Area  vegetation  includes  grasses,  shrubs,  lichens, 
and  mosses. 

The  landform  at  Imnaviat  Creek  is  mostly  a  hill  slope 
(more  than  20%  slope),  with  a  portion  of  the  area  in  a 
basin.  The  aspect  of  the  slope  ranges  from  west  to  north. 
Soils  of  Imnaviat  Creek  consist  of  about  10  cm  of  live 
and  dead  organic  material  over  5  to  10  cm  of  partially 
decomposed  organic  material  mixed  with  silt,  all  over 
glacial  till,  and  are  listed  as  Histic  Pergelic  Cryaquepts 
(Hinzman  et  al.  1991).  Vegetation  in  Imnaviat  consists 
of  mostly  water-tolerant  plants  such  as  tussock-sedges, 
mosses  associated  with  lichens,  and  shrubs  such  as 
dwarf  birch,  alder,  and  willows  (Kane  et  al.  1994).  Tus¬ 
sock-sedges  are  common  on  the  ridge  tops  and  on 
slopes,  but  lacking  in  the  valley  bottom.  Mosses  are 
dominant  in  the  valley  bottom  where  the  environment 
is  wettest. 

The  Toolik  Lake  site  is  contained  on  a  hillslope  (up 
to  20%  slope  with  an  aspect  ranging  from  east  to  north) 
on  colluvial  deposits  and  a  portion  of  the  area  is  on 


frozen  water  (pond).  More  than  half  of  the  area  con¬ 
sists  of  moist  acidic  tundra  (soil  pH  <  5.0)  dominated 
by  tussock-sedges,  dwarf  shrubs,  and  mosses.  The 
remaining  area  is  moist  non-acidic  (soil  pH  >  5.0)  and 
dry  acidic  tundra  vegetation  (Walker  1998).  The  soils 
in  Toolik  Lake  are  classified  as  a  coarse-loamy,  mixed 
Pergelic  Cryochrepts  composed  of  layers  of  peat,  loam, 
and  sandy  loam  (Ping  et  al.  1994). 

West  Dock  is  at  an  elevation  of  3  m  above  mean  sea 
level.  The  soil  is  composed  of  layers  of  peat,  loamy 
sand,  and  fine  sandy  loam  (Ping  et  al.  1994)  and  is  classi¬ 
fied  as  coarse-loamy,  mixed  Pergelic  Ruptic-Histic 
Cryaquept.  The  vegetation  consists  of  sedge  shrubs, 
mosses,  and  lichens. 

Betty  Pingo  and  West  Dock  are  located  in  a  wetland 
complex  with  little  topographic  relief  near  Prudhoe  Bay 
and  are  considered  part  of  the  upper  Kuparuk  River 
drainage  basin  (Walker  et  al.  1980).  West  Dock,  Betty 
Pingo,  Happy  Valley,  Imnaviat  Creek,  and  Toolik  Lake 
are  located  on  the  Kuparuk  River  drainage  basin  (Nel¬ 
son  et  al.  1997). 

Fort  McCoy,  Wisconsin 

In  contrast  to  the  Alaska  sites,  this  site  is  in  a  non- 
permafrost  area.  The  average  freezing  season  is  1 1 7  days 
and  expected  frost  depth  is  0.50  m.  Thaw  depths  were 
collected  on  23  March  1995,  at  Fort  McCoy,  Wiscon- 
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Figure  3.  Location  map  for  Minnesota  and  Wisconsin  sampling  sites. 
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Figure  4.  Sampling  grid  for  Fort  McCoy. 


sin  (44°N,  90°W),  shown  in  Figure  3,  along  a  narrow 
580-m-long  trail  (Kestler  1996,  Kestler  et  al.  1999).  The 
trail  consisted  of  a  saturated,  thawing  soil  layer  over 
frozen  ground.  Thaw  depths  were  measured  at  7.6-m 
intervals  along  the  trail  centerline  and  at  30.50- 
m  intervals  in  the  right  and  left  wheel  paths,  which 
were  approximately  1.2  m  off  the  trail  centerline. 
Figure  4  is  the  typical  sampling  pattern.  A  metal 
rod  was  used  to  determine  thaw  depth  by  push¬ 
ing  the  rod  into  the  ground  until  it  met  resistance 
and  by  checking  the  temperature  at  this  interface. 

The  soil  on  this  trail  was  classified  as  silty 
sand.  The  trail  water  content  ranged  from  8  to 
34%  with  an  average  of  18.6%  and  standard  de¬ 
viations  of  4.4.  The  optimum  moisture  content 
(moisture  of  the  soil  when  compacted  at  maxi¬ 
mum  density)  of  the  soil  mineral  ranged  from  8.2 
to  1 1.4%.  The  trail  consisted  of  a  thin  (0-  to  10- 
cm)  vegetated  layer  on  top  of  the  mineral  soil. 

The  percent  slope  of  the  terrain  ranged  from  flat 
to  gently  sloping  (with  a  maximum  of  one  short 
section  of  10%  slope).  Trees  were  present  along 
the  sides  of  the  trail,  contributing  to  variations  in 
shading  on  the  trail  surface. 

Swan  Lake,  Minnesota 

Thaw  depths  were  collected  on  a  shallow  de¬ 
pression  in  the  Swan  Lake  watershed  (Fig.  3), 
near  Morris,  Minnesota  (45°41'N,  95°48'W),  by 


Dr.  Brenton  Sharratt  of  the  U.S.  Department  of  Agri¬ 
culture  (USDA),  North  Central  Soil  Conservation  Re¬ 
search  Laboratory  (Sharratt  et  al.  1999).  The  elliptical 
depression  or  prairie  pothole  was  about  2  ha  in  area 


Figure  5.  Sampling  grid,  site  contour,  and  soil  map  for  Swan 
Lake  Watershed. 
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Figure  6.  FERF  grid  sampling  patterns. 


with  2  m  of  relief.  The  major  axis  of  the  depression  is 
oriented  NW-SE.  Thaw  depth  measurements  were  taken 
along  both  the  major  and  minor  axis  (Fig.  5)  using  frost 
tubes  during  spring  1991  and  1992. 

The  soil  types  present  within  the  depression  included 
loam,  silty  clay  loam,  and  clay  loam.  Wheat  was  grown 
on  this  site  in  the  spring  of  1990  and  1991.  The  depres¬ 
sion  was  chisel-plowed  after  fall  harvest  (October)  each 
year.  Unlike  the  Alaska  sites,  this  site  is  without  perma¬ 
frost.  The  maximum  soil  frost  depth  of  this  site  is  1 .2 
m,  and  average  freezing  period  is  130  days. 

FERF 

Thaw  depths  were  collected  in  a  controlled  and 
prepared  freeze-thaw  experiment  in  the  FERF  on  a 
homogeneous  soil  compacted  at  specified  moisture 
content  and  density.  This  site  is  an  indoor,  temperature- 
controlled  environment  that  should  reflect  homogen¬ 
ous  soil  conditions  and  eliminates  the  climatic,  terrain, 
and  solar  radiation  effects  of  the  natural  environment. 
Thaw  depths  were  measured  over  the  100-square-meter 
area  by  pushing  a  metal  rod  into  the  soil  and  were  veri¬ 
fied  by  digging  to  the  frost  layer.  Figure  6  shows  typi¬ 
cal  grid  spacings  used  for  these  measurements.  (Some 
of  the  data  sets  were  collected  on  a  tighter  grid  than 
others.) 

The  soil  type  was  characterized  as  silt.  The  dry  den¬ 
sity  of  the  soil  varied  from  1.39  to  1.70  g/cm3  and  the 
moisture  content  ranged  from  20  to  30%  (optimum 
water  content  is  22%).  The  soil’s  maximum  density 
from  the  standard  proctor  test  was  1.52  g/cm3  (Shoop 
et  al.  2000). 

STATISTICAL  ANALYSIS 

Thaw  depths  were  analyzed  to  assess  the  distribu¬ 
tion  of  thaw  as  well  as  to  determine  the  spatial  effects. 
The  statistical  calculations  include  the  average,  stan¬ 


dard  deviation  and  coefficient  of  variation  (CV).  The 
ratio  of  standard  deviation  to  the  mean,  eq  1 ,  is  often 
used  to  describe  the  distribution  (i.e.,  the  histogram) 
(Isaaks  and  Srivastava  1989).  A  CV  greater  than  one 
indicates  the  presence  of  some  erratic  sample  values 
and  high  variability  in  the  data.  A  CV  less  than  one 
indicates  low  variability  in  the  data  set.  Histograms  of 
the  thaw  depth  are  displayed  for  each  of  the  data  sets. 

CV  =  -  CD 

m 

where  CV  =  coefficient  of  variation 
o  =  standard  deviation 
m  =  mean. 

Alaska  North  Slope 

Representative  histograms  from  each  of  the  sites  in 
Alaska  are  shown  in  Figure  7  (all  histograms  are  given 
in  Appendix  A).  Statistical  analysis  for  each  of  the  thaw 
depth  measurement  events  for  all  data  sets  is  given  in 
Tables  2  and  3.  As  mentioned  earlier,  there  were  miss¬ 
ing  thaw  depths  on  the  grid  because  the  point(s)  or  por¬ 
tion  of  the  grid  was  located  on  lakes,  standing  water, 
rocks,  and  roads.  Thus  code  200  was  used  for  lakes  and 


Table  2.  Thaw  depth 

coefficient-of-variation 

distribution  for  all  sites. 

Sites 

CV  Range 

Average  CV 

Atquasuk,  Alaska 

0.4-0. 5 

0.43 

Barrow,  Alaska 

0.3-0. 5 

0.32 

Betty  Pingo,  Alaska 

0.3 

0.33 

Happy  Valley,  Alaska 

0. 2-0.4 

0.28 

Imnaviat  Creek,  Alaska 

0.2-0. 3 

0.25 

Toolik  Lake,  Alaska 

0. 3-0.4 

0.32 

West  Dock,  Alaska 

0.3 

0.26 

Fort  McCoy,  Wisconsin 

0.5 

0.45 

Swan  Lake,  Minnesota 

0.2-0. 6 

0.39 

FERF,  Hanover,  N.H. 

0. 1-0.6 

0.21 
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Figure  7.  North  Slope  thaw-depth  histogram. 


standing  water  and  300  was  used  for  rocks,  gravel, 
roads,  etc.  These  codes  were  culled  for  our  analysis. 
All  seven  sites  have  a  very  low  coefficient  of  variation, 
ranging  from  0.2  to  0.5,  indicating  minimum  variation 
in  the  data  sets. 

The  histograms  for  Atquasuk,  Betty  Pingo,  and  West 
Dock  are  somewhat  skewed  with  a  long  tail  of  high 
thaw-depth  values.  The  skewed  data  sets  were  located 
in  the  coastal  areas  with  relatively  flat  terrain;  howev¬ 
er,  it  is  unclear  how  these  factors  may  have  affected  the 
distribution.  Thaw  depth  for  Barrow,  Happy  Valley, 
Imnaviat  Creek,  and  Toolik  Lake  more  closely  resem¬ 
ble  a  normal  distribution.  These  behaviors  are  consis¬ 
tent  for  all  data  sets  for  each  site  (App.  A). 


Fort  McCoy,  Wisconsin 

A  histogram  for  thaw  depths  on  the  trail  at  Fort 
McCoy  showed  scattered  distribution  with  a  wide  range 
and  CV  of  0.5  as  shown  in  Figure  8.  Table  3  shows  the 
minimum  and  maximum  thaw  depth.  As  the  data  were 
collected  for  only  one  date,  no  temporal  comparisons 
are  possible. 

Swan  Lake,  Minnesota 

The  coefficient  of  variation  for  the  13  data  sets  from 
Swan  Lake,  Minnesota,  varied  from  0.2  to  0.6  (Tables 
2  and  3).  Figure  9  illustrates  typical  histograms  for  the 
data  sets  and  all  of  the  data  histograms  are  shown  in 
A8.  The  thaw  depths  showed  spread  distribution.  These 


7 


to  contents 


0.00  20.00  40.00  60.00  80.00 

Thaw  Depth  (cm) 

Figure  8.  Fort  McCoy  thaw-depth  histogram. 


thaw  depth  variations  were  due  to  the  in¬ 
fluence  of  variations  of  solar  radiation  and 
topography  (slope  and  aspect).  Because 
of  the  shallow  depression,  the  southern- 
facing  slope  gains  more  intense  sunlight 
during  part  of  the  day  while  the  northern- 
facing  slope  gets  the  self-shadowing  ef¬ 
fect.  Sharrat  et  al.  (1999)  mentioned  that 
rapid  thawing  at  location  9  and  possibly 
location  5  occurred  when  water  was 
present  in  the  center  of  the  depression.  The 
progression  of  thaw  with  time  is  shown 
in  Figure  10,  showing  anomalous  thaws 
at  locations  5,  9,  and  20.  In  general  the 
thaw  depth  progressed  downward  evenly 
with  time. 

FERF 

The  average  thaw  depth  for  the  FERF 
data  sets  ranged  from  2.5  to  20.8  cm.  Stan¬ 
dard  deviations  ranged  from  1  to  3.5  cm 
(Tables  2  and  3).  For  all  data  sets,  the  thaw 
depth  is  higher  at  one  end  of  the  test  basin 
(Fig.  11).  This  is  true  for  most  of  the  data 
sets  in  the  FERF,  and  is  probably  due  to 
changes  in  the  air  temperature  and  frost 
properties  at  the  north  end  of  the  test  sec¬ 
tion  where  the  floor  of  the  test  basin  ramps 
up  toward  the  building  entrance.  Figure 
12  shows  the  histograms  for  three  of  the 
thaw  events. 

COMPARISON  OF  SITES 


Thaw  Depth  (cm) 


A  comparison  of  variability  at  each  site  Figure  9.  Typical  histograms,  Swan  Lake. 

is  given  in  Table  2.  The  FERF  data  sets 
showed  the  lowest  CV,  with  an  average 
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Figure  10.  Elevation  versus  thaw  depth  at  Swan  Lake. 


9 


to  contents 


-2/5/91 

-2/6/91 

2/8/91 

3/18/91 

-3/20/91 

-3/21/91 

-3/22/91 

-3/25/91 

3/26/91 


Elevation  NE  Aspect 


Figure  10  (cont’d).  Elevation  versus  thaw  depth  at  Swan  Lake. 
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Figure  12.  Typical  histograms,  FERF. 


CV  of  0.2 1 .  This  value  can  be  considered  indicative  of 
a  bare  (no  vegetation),  homogenous  soil  in  a  tempera¬ 
ture-controlled  environment,  with  no  solar  input.  The 
CV  for  Fort  McCoy,  Wisconsin,  was  the  highest  due  to 
trees  present  along  the  sides  of  the  trail  (which  contrib¬ 
uted  a  disparity  in  shading)  and  to  water  ponding  and  a 
stream  crossing  in  some  parts  of  the  trail.  The  variation 
of  thaw  depth  distribution  on  Swan  Lake,  Minnesota, 
data  sets  is  due  to  the  slope  and  aspect  of  the  site  (Shar- 
rat  et  al.  1999).  Atquasuk  has  a  higher  CV  than  the  rest 
of  the  Alaska  sites  but  the  reasons  are  unclear. 

SPATIAL  STATISTICS 

Variogram  functions  to  determine  the  correlation 
distance  for  thaw  depth  were  used  to  analyze  the  data 
sets.  The  variogram  function  is  a  mathematical  model 


that  expresses  the  statistical  variance  as  a  function  of 
the  separation  distance.  For  example,  several  thaw  depth 
measurements  taken  at  nearly  the  same  point  (e.g.,  point 
20)  will  have  nearly  the  same  value,  so  the  variation 
will  be  small.  If  spatial  correlation  exists,  sequentially 
increasing  the  separation  distance  between  data  points 
will  yield  a  larger  variance  in  the  data.  This  is  the  pro¬ 
cess  used  to  generate  a  variogram  where  separation  dis¬ 
tance  is  plotted  on  the  x-axis  and  the  variance  on  the 
y-axis.  Eq  2  defines  the  variogram.  If  the  variance 
increases  with  the  separation  distance,  then  the  thaw  depth 
is  correlated  with  distance.  At  the  point  when  an  increase 
in  separation  distance  no  longer  causes  a  correspond¬ 
ing  increase  in  the  variance,  the  variogram  reaches  a 
plateau.  The  distance  corresponding  to  the  start  of  the 
plateau  is  called  the  range  and  the  variance  value  at  the 
plateau  is  called  the  sill.  The  vertical  jump  of  variance 
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Figure  14.  Three  most  commonly  used  transition  models  with  the  same  range 
and  sill  (spherical,  exponential,  and  Gaussian). 


Figure  15.  Power  model. 
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at  the  origin  when  separation  distance  is  zero  is  called 
the  nugget  effect.  Measurement  error  and/or  natural 
occurrences  can  cause  this  discontinuity  at  the  origin. 
These  values  are  illustrated  in  Figure  13. 
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The  basic  variogram  models  are  divided  into  two 
types;  those  that  reach  a  plateau,  which  are  called  tran¬ 
sition  models,  and  those  that  do  not.  Theoretical  models 
such  as  spherical,  exponential,  Gaussian,  and  power  (or 
linear)  models  are  used  to  fit  the  thaw  depth  on  the  vario¬ 
gram  plot  (Fig.  14  and  15).  The  spherical  model  is  the 
most  commonly  used  transition  model  (eq  3)  and  has  a 
linear  behavior  at  small  separation  distance  near  the 
origin  but  flattens  out  at  larger  distance  and  reaches  a 
sill  at  a.  The  tangent  at  the  origin  reaches  the  sill  at 
about  two-thirds  of  the  correlation  range.  The  exponen¬ 
tial  model  is  another  commonly  used  transition  model 
(eq  4).  The  exponential  model  rises  more  steeply  near 
the  origin  and  flattens  out  more  gradually.  Like  the  expo¬ 
nential  model,  the  Gaussian  model  reaches  its  sill 
asymptotically  with  a  parabolic  behavior  near  the  ori¬ 
gin,  and  the  sill  is  defined  as  the  “practical  range”  or 
the  distance  at  which  the  variogram  value  is  95%  of  the 
sill  (eq  5).  The  power  or  linear  model  is  not  a  transition 
model  since  it  does  not  reach  a  sill  but  continuously 
increases  with  separation  distance  h  (eq  6).  The  soft¬ 
ware  we  used  for  variogram  analysis  is  Variowin  (Pan- 


natier  1 996).  According  to  ASTM  standards,  variogram 
analysis  should  have  at  least  20  pairs  of  data  (ASTM  D 
1996). 


Spherical  Variogram  Model 
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Gaussian  Variogram  Model 
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Alaska  North  Slope 

The  variogram  model  that  best  fit  the  North  Slope 
thaw  depth  data  was  the  spherical  model.  Spatial  cor¬ 
relations  were  found  on  Atquasuk,  Betty  Pingo,  Happy 
Valley,  and  West  Dock  data  sets.  These  data  grids  are 
situated  on  the  coastal  plain  except  for  Happy  Valley 
(near  the  foothills).  Atquasuk  shows  the  correlation 
range  of  approximately  660  m  and  with  sill  ranging 
from  230  to  440  depending  on  the  sampling  date  (Fig. 
1 6a).  The  correlation  range  for  the  three  sampling  dates 
is  constant.  Betty  Pingo  exhibits  a  correlation  range  of 
up  to  450  m  and  sill  varying  from  140  to  370  (Fig.  16c). 
The  thaw  depth  correlation  range  for  Betty  Pingo  data 


a.  Atquasuk.  b.  Barrow. 

Figure  16.  Variogram  plots  for  Alaska  North  Slope. 


13 


to  contents 


c.  Betty  Pingo.  d.  Happy  Valley. 


e.  Imnaviat  Creek. 


f.  Toolik  Lake. 


h 


g.  West  Dock. 


Figures  16  (cont’d).  Variogram  plots  for  Alaska  North  Slope. 
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Figure  17.  Fort  McCoy  variogram. 


Figure  18.  Swan  Lake  variogram. 


sets  increases  for  the  first  two  sampling  dates  and 
remains  constant  for  the  third  sampling  date.  The  first 
thaw  depth  measurements  (8  July)  on  Happy  Valley 
showed  no  correlation;  however,  the  following  two 
measurements  showed  a  correlation  range  of  250  and 
225  m  with  sill  of  110  and  82,  respectively  (Fig.  16d). 
The  sill  for  West  Dock  ranged  from  1 10  to  185  with  a 
nearly  constant  range  of  400  m  (Fig.  16g). 

The  variogram  plots  for  Barrow  and  Toolik  Lake 
data  sets  displayed  correlation  for  all  data  points.  The 
sill  and  correlation  range  values  were  undetermined. 
Barrow  exhibits  a  linear  increase  in  variance  with  sep¬ 
aration  distance  (Fig.  16b).  For  Barrow  a  much  larger 
sampling  area  is  needed  to  discern  correlation  at  these 
separation  distances.  Toolik  Lake  variogram  also 
showed  a  gradual  increase  in  variance  starting  at  100 
m  with  a  range  that  cannot  be  determined  (Fig.  16f); 
the  farther  apart  the  points,  the  greater  the  variance. 
Similar  to  Barrow,  a  larger  sampling  area  is  needed  for 
Toolik  Lake  to  detect  the  spatial  correlation  range.  Both 
Toolik  Lake  and  Happy  Valley  are  near  each  other  and 
situated  on  the  foothills,  which  have  similar  vegetation, 
slopes,  and  local  weather. 

Figure  16e  shows  the  Imnaviat  Creek  variogram 
where  the  statistical  variance  fluctuated  throughout  the 
sampling  area  with  no  correlation.  The  total  absence  of 
spatial  correlation  for  Imnaviat  Creek  meant  that  the 
separation  distance  or  sampling  interval  of  100  m  is  so 
far  apart  that  one  sampling  location  is  unrelated  to  the 
next  at  these  sites.  It  also  could  be  due  to  the  terrain 
slope,  vegetation,  or  local  weather. 

Fort  McCoy,  Wisconsin 

The  thaw  depths  that  were  measured  at  Fort  McCoy 
had  a  correlation  range  of  6 1  m  and  a  maximum  statis¬ 
tical  variance  (sill)  of  22  (Fig.  17).  The  spherical  vari¬ 
ogram  model  is  the  best  fit  for  these  data.  This  residt  is 
consistent  with  a  similar  analysis  by  Kestler  (1996). 


Swan  Lake,  Minnesota 

The  variogram  plots  for  Swan  Lake  did  not  reveal 
any  spatial  correlation.  Figure  18  shows  a  typical  vari¬ 
ogram  plot  with  scattered  points.  The  layout  (Fig.  5) 
with  a  sampling  distance  of  15  m  is  too  far  apart,  so 
that  correlation  distance  was  undetectable.  Perhaps  a 
shorter  sampling  distance  or  a  regular  sampling  grid  is 
needed  to  determine  the  correlation.  All  the  variogram 
plots  for  Swan  Lake  data  sets  are  shown  in  Appendix 
B 1 .  Solar  radiation  accounted  for  1 1  %  of  the  variabili¬ 
ty  in  thaw  depth  across  the  depression,  with  more  rapid 
thaw  occurring  at  locations  with  greater  radiation 
because  of  the  geometry  of  the  terrain  (slope  and  aspect). 
Most  of  the  lower  elevation  thawed  faster  than  the  high¬ 
er  areas  (Sharratt  et  al.  1999). 

FERF 

Best  fits  were  obtained  using  exponential,  power, 
and  linear  variogram  models  for  FERF  data  sets.  The 
typical  results  from  variogram  analysis  are  shown  in 
Figure  19a  and  b.  Plots  for  all  the  data  sets  are  located 
in  Appendix  B2.  All  plots  showed  that  the  thaw-depth 
measurements  have  an  increase  in  variance  with  separa¬ 
tion  distance;  however,  the  correlation  range  cannot  be 
determined.  Because  the  distribution  of  the  thaw  depth 
gradually  slopes  across  the  test  section  as  shown  in  Fig¬ 
ure  1 1 ,  the  variogram  plots  tend  to  have  a  concave  corre¬ 
lation.*  The  data  sets  with  the  largest  number  of  sam¬ 
pling  points  most  strongly  show  this  concave  behavior. 
The  absence  of  the  sill  indicates  that  the  correlation 
range  was  larger  than  the  maximum  spatial  distance  of 
27  m  and  a  larger  sampling  area  is  needed  to  determine 
the  correlation  distance. 


*  Personal  communication,  Dr.  Matthew  Sturm,  CRREL,  Fairbanks, 
Alaska,  2000. 


15 

to  contents 


a.  21  June  1992.  b.  19  December  1991. 

Figure  19.  Typical  variogram  plots  of  FERF  data  sets. 


APPLICATIONS 

Variogram  analysis  allows  one  to  quantify  the  spatial 
variability  of  the  thaw  depth  for  various  types  of  ter¬ 
rain.  Standard  statistics  (i.e.,  standard  deviation,  histo¬ 
gram)  do  not  relate  to  the  relative  positions  of  the  test 
data  but  are  essential  information  on  the  variation  of 
thaw  depth.  The  correlation  range  and  variation  of  thaw 
depth  over  a  terrain  can  be  used  to  quantify  the  thaw 
depth  variability  in  models  where  thaw  depth  is  needed, 
such  as  mobility  performance  models. 

The  depth  of  thaw  causes  a  significant  impact  on 
vehicle  mobility  (Shoop  1990).  Soils  with  a  thawed 
layer  can  be  saturated  due  to  snowmelt,  rain,  and  water 
drawn  to  the  soil  during  the  freezing  process.  The  under¬ 
lying  frozen  layer  traps  this  moisture,  so  the  upper  soil 
layer  is  wet  and  weak.  A  thaw-depth  algorithm  has  been 
incorporated  into  a  vehicle  performance  model  called 
NATO  Reference  Mobility  Model  II  (NRMMII)  (Rich¬ 
mond  et  al.  1995).  NRMMII  predicts  maximum  attain¬ 
able  speed  of  a  vehicle  operating  on  a  terrain  (Ahlvin  and 
Haley  1992).  The  terrain  is  divided  into  a  terrain  unit 
for  corresponding  soil  type,  soil  strength,  soil  wetness 
index,  soil  strength,  obstacles,  and  surface  roughness. 
Thaw  depth  is  an  input  variable  to  a  terrain  or  scenario 
file  for  NRMMII  and  is  assumed  uniform  throughout 
the  terrain  or  can  be  specified  for  each  terrain  unit. 

Kestler  (1996)  used  statistical,  geostatical  (vario¬ 
gram),  and  regression  analyses  to  establish  whether 
variation  in  subgrade  reflects  the  performance  of  a  test 
section.  Her  study  concluded  that  thaw  depth,  water 
content,  percentage  of  water  standing,  and  soil  strength 
influenced  the  test  section  performance.  Also,  Kestler 
et  al.  (1994)  identified  spatial  correlation  as  a  function 
of  distance  using  the  falling  weight  deflectometer 
(FWD)  data  to  assess  the  overall  pavement  structure. 
The  FWD  is  towed  by  a  vehicle  and  is  used  to  deter¬ 


mine  pavement  stiffness  and  deflection.  Using  the 
knowledge  of  the  lengths  over  which  any  particular 
indicator  value  is  correlated  can  translate  into  substan¬ 
tial  savings  by  optimizing  pavement  evaluation  tests. 

Temporal  and  spatial  variability  on  soil  strength  has 
been  investigated  for  off-road  mobility  to  define  the 
soil  strength  uncertainties,  and  for  forecasting  soil 
strength  (Mason  et  al.  1997).  Mason  et  al.  used  soil 
strength  data  to  develop  coefficients  to  forecast  the 
changes  of  soil  strength  with  time,  and  spatial  relation¬ 
ship  on  soil  strength  with  the  terrain  elevation  using 
variogram  and  kriging  methods. 

DISCUSSION 

The  average  thaw  depth  versus  the  standard  devia¬ 
tion  of  all  the  data  sets  from  all  sites  is  shown  in  Figure 
20.  The  variance  in  the  thaw  is  proportional  to  the  thaw 
depth  according  to  eq  7  and  8. 


°thaw  0.3007  WJthaw 

(7) 

°thaw  —  0.1513  W?thaw 

(8) 

where  Othaw  is  the  standard  deviation  of  the  thaw  depth 
and  »?thaw  is  the  average  thaw  depth  (cm).  The  controlled 
test  section  in  the  FERF  has  a  slightly  smaller  slope. 

At  sites  where  data  was  collected  with  thaw  pro¬ 
gression  (i.e.,  Alaska),  spatial  correlation  ranges  are 
fairly  constant  as  the  thaw  depth  progresses.  In  Swan 
Lake,  Minnesota,  data  was  collected  with  time  but  no 
correlation  range  was  found  as  discussed  previously. 

The  data  sets  were  collected  at  various  sampling 
scales.  For  example,  at  Swan  Lake,  the  data  were  mea¬ 
sured  within  the  confines  of  the  shallow  depression, 
and  in  the  case  of  the  FERF,  it  is  the  dimension  of  the 
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Figure  20.  Variability  with  thaw  depth. 


installation.  These  limit  the  spatial  variability  analysis 
to  determine  the  correlation  range  on  these  sites,  and  are 
apparent  in  the  variograms. 

Solar  radiation  is  one  of  the  influence  factors  for 
slopes  with  varying  orientation  and  gradient.  Thus,  the 
change  in  solar  radiation  accounts  for  some  of  the  varia¬ 
bility  in  the  Alaska  sites  on  the  foothills  (Nelson  et  al. 
1997),  and  in  the  Swan  Lake,  Minnesota,  site  (Sharrat  et 
al.  1999). 

The  correlation  range  on  Table  2  can  be  used  as  an  ini¬ 
tial  approximation  for  the  specific  terrain.  For  example,  in 
the  case  for  seasonal  thaw  depth  without  permafrost,  the 
correlation  range  can  be  estimated  at  about  6 1  m. 

These  results  can  be  used  for  analysis  of  the  impact 
of  thaw  distribution  on  predictive  models,  such  as  for 
predicting  vehicle  mobility.  Future  research  should 
include  evaluation  of  impact  of  climate,  vegetation,  sur¬ 
face  roughness,  solar  radiation,  and  slope  on  thaw  vari¬ 
ability,  and  how  they  change  with  time. 

SUMMARY 

Statistical  and  spatial  analyses  were  performed  on 
existing  thaw  depth  data  sets  from  a  variety  of  climatic 
and  terrain  conditions  to  determine  the  variability  of 
thaw.  These  data  include  permafrost  sites  in  Alaska,  tem¬ 
perate  sites  in  the  continental  United  States,  including 
sites  in  Minnesota  and  Wisconsin,  and  a  large-scale, 


controlled  test  section  with  no  solar  input  in  the  Frost 
Effects  Research  Facility  (FERF)  at  the  Cold  Regions 
Research  and  Engineering  Laboratory  (CRREL)  in  Flan- 
over,  New  Flampshire. 

Variation  in  elevation,  slope,  aspect,  vegetation,  and 
solar  radiation  contribute  to  the  variation  in  thaw.  Coef¬ 
ficients  of  variation  ranged  from  very  low  values  of  0. 1 
at  the  controlled  site  (FERF)  to  0.5  and  0.6  at  sites  with 
non-uniform  elevation,  slope/aspect,  vegetation,  or  sur¬ 
face  water.  Spatial  statistic  analysis  of  the  thaw  depth 
data  show  that  the  correlation  range  at  each  site  stays 
constant  as  thaw  progresses.  The  correlation  range  at 
each  site  can  be  used  as  a  first  approximation  for  esti¬ 
mating  correlation  of  similar  sites.  The  terrain  and  sur¬ 
face  energy  factors  were  not  specifically  addressed, 
making  broad  generalization  regarding  the  impacts  of 
specific  factors  difficult. 

Even  with  the  wide  differences  among  sites,  the  stan¬ 
dard  deviation  tends  to  increase  smoothly  with  thaw 
depth.  A  compilation  of  all  data  (Fig.  20)  shows  an 
increase  according  to 

°thaw  ~  0.03007  WJthaw-  (9) 

Most  data  fit  this  trend  with  the  controlled  site  having  a 
slightly  smaller  slope.  Because  the  Alaska  sites  were 
underlain  by  permafrost,  most  of  the  deepest  thaw 
depths  were  collected  in  this  environment. 
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Results  can  be  used  to  spatially  distribute  soil  prop¬ 
erties  based  on  point  data  or  one-dimensional  models, 
or  to  populate  sparse  data  sets  with  terrain  properties. 
They  are  also  useful  for  analysis  of  impact  of  thaw  dis¬ 
tribution  on  predictive  models,  such  as  for  predicting 
vehicle  mobility. 
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Figure  A1 .  Atquasuk,  Alaska. 


Figure  A2.  Barrow,  Alaska. 
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Figure  A3.  Betty  Pingo,  Alaska. 
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Figure  A4.  Happy  Valley,  Alaska. 
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Figure  A5.  Imnaviat  Creek,  Alaska. 
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Figure  A7.  West  Dock,  Alaska. 
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Figure  A8.  Swan  Lake,  Minnesota. 
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Figure  A8  (cont’d). 
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Figure  A9.  FERF 
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APPENDIX  B1.  SWAN  LAKE  THAW-DEPTH  VARIOGRAMS 
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